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Abstract 
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o\ 

O^ j We present measurements of the production cross sections times leptonic 

W) branching fractions and the transverse momentum distributions of W and Z 

bosons in pp collisions at = 1.8 TeV using data collected with the D0 
detector at the Fermilab Tevatron pp collider. A preliminary measurement of 

£S) the W charge asymmetry is also presented. 
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At y/s = 1.8 TeV, production of W and Z bosons in pp collisions proceeds primarily via qq 
annihilation accompanied by an initial state gluon radiation which produces the transverse 
momentum, p?, of W and Z bosons. Absolute predictions for the inclusive production cross 
sections, aw and az, have been calculated to order a 2 s by van Neerven et al. [1]. In the 
low px region (px ,Px < 20 GeV/c) multiple soft gluon emission is expected to dominate 
the initial state radiation and the differential cross section, da/dpx } has been calculated 
using a soft gluon resummation technique [2-6]. In the high pj region (px ,Px > 20 GeV/c) 
perturbative QCD calculations are expected to describe the da / dpx [7]. Thus measurements 
of the inclusive and differential cross sections of W and Z production provide tests of QCD. 
Experimentally, use of leptonic decays of W and Z bosons, which do not involve final state 
strong interactions, allows for high precision measurements of their inclusive processes. 

In this report we present measurements of the production cross sections times leptonic 
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branching fractions and the px distributions of W and Z bosons using data collected with the 
D0 detector [8] in the 1992-1993 run at the Fermilab Tevatron pp collider at y/s = 1.8 TeV. 
We present preliminary results using a partial data sample from the 1994-1995 run. We 
also present a preliminary measurement of the W charge asymmetry using W — > \iv sample, 
from which the information on parton distribution function (pdf) can be extracted. 



where N f, s is the number of observed events, N^kgd is the number of expected background 
events, Aw is the kinematic and geometric acceptance, ew is the detection efficiency, and C 
is the integrated luminosity used in the analysis. The Z boson cross section, oz • B(Z — > //), 
is calculated in a similar fashion. 

Electrons were detected in hermetic, uranium liquid-argon calorimeters with an energy 



resolution of about 15%/y -E'(GeV). The central and end calorimeter regions were used in 
both the W and Z analyses, covering pseudorapidity (77) range: \rj\ < 1.1 and 1.5 < \r)\ < 2.5, 
respectively. Muons were detected as tracks in three layers of proportional drift tube cham- 
bers outside the calorimeter: one 4-plane layer is located inside a magnetized iron toroid 
and two 3-plane layers are located outside, providing coverage for \rj\ < 3.3. The muon mo- 
mentum resolution is ail/p) = 0.18(p — 2)/p 2 © 0.008 (with p in GeV/c). Muons that 
passed through the central iron toroid (\rj\ < 1.0) were used in the cross section measure- 
ments. Neutrinos were inferred from the observed missing transverse energy (Et) which was 
calculated using all the energy detected in the calorimeter cells out to pseudorapidity of 4.2. 
For the electron channel decays, the Et resolution was dominated by the underlying event 
and is ~ 3 GeV. For the muon channel decays, the muon transverse momentum was added 
to the calorimeter energy to calculate the total Et } and the muon momentum resolution 
dominated the Et resolution. 



THE INCLUSIVE PRODUCTION CROSS SECTIONS [9] 



The W boson inclusive cross section is calculated as 
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TABLE I. Estimates of Backgrounds 
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1.2 ± 0.1 




1.7 ± 0.3 


Total Background(%) 


5.7 ± 0.5 


4.0 ± 1.4 


22.1 ± 1.9 


10.1 ± 3.7 



The W and Z electron channel analyses based their event selection on a sample ob- 
tained with a single electron trigger (Et > 20 GeV). Offline, it was required that there be 
at least one electron with Et > 25 GeV that passed "tight" electron identification cuts. 
Details of the electron identification are given in Ref. [10], with the main features being 
an electromagnetic (EM) cluster in the calorimeter with a matching track in the central 
tracking chambers. The electron was required to have the isolation variable less than 0.1, 
where the isolation is defined as 7=(£ , tot (0.4)-£ , EM(0.2))/£ , EM(0.2), and £' tot (0.4) is the total 
calorimeter energy inside a cone of radius \/ Ar] 2 + Ac/) 2 = 0.4 and £^em(0.2) is the elec- 
tromagnetic energy inside a cone of 0.2. The cluster was also required to have transverse 
and longitudinal shapes consistent with those expected for an electron based on test beam 
measurements and Monte Carlo simulations. To select W — > ev candidates, in addition to 
the "tight" electron with Et > 25 GeV, events were required to have missing transverse 
energy Et > 25 GeV. To select Z — > ee candidates, in addition to the "tight" electron with 
Et > 25 GeV, events were required to have a second electron with Et > 25 GeV but the 
electron identification requirements were loosened by not requiring the track match in order 
to increase the efficiency. The invariant mass of the electron pair was required to be in the 
range 75 < M ee < 105 GeV/c 2 . In an analysis of the 1992-1993 data sample, corresponding 
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FIG. 1. Mass spectra from the 1992-1993 run. The points are the data, the shaded areas are 
the estimated backgrounds, and the histograms are the sum of the MC predictions and estimated 
backgrounds. 

to 12.8 ± 0.7 pb -1 , we found 10338 W and 775 Z candidate events. The mass spectra for 
the W — > ev and Z — > ee events are shown in Fig. 1. 

The muon channel W and Z analyses used an event sample obtained with a single muon 
trigger (px > 15 GeV). Offline, the events were required to have a reconstructed muon 
with px > 20 GeV. For W — > \iv events, the missing transverse energy was required to be 
Jjjf > 20 GeV. For Z — > fi/j, events, the offline threshold on the second muon was lowered 
to 15 GeV and the muon identification criteria were loosened. The main features of the 
muon identification (see Refs. [9,10] for details) include a good quality muon track that 
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FIG. 2. Mass spectra (electron decay channels above, muon decay channels below) from a 
partial sample of the 1994-1995 data. 

has a calorimeter confirmation signal and has a stringent match with a track in the central 
detector. Cosmic ray background was reduced by rejecting muons that also had hits or tracks 
within 10° in 6 and 20° in <f> in the muon chambers on the opposite side of the interaction 
point. For the W — > \iv selection, events that were Z — > fi/j, candidates were removed. From 
an analysis of the 1992-1993 data sample, corresponding to 11.4±0.6 pb -1 , we found 1665 W 
and 77 Z candidate events. The observed mass spectra for the W — > \iv and Z — > fi/j, events 
are shown in Fig. 1. 

A preliminary analysis of a partial sample of the 1994-1995 data, corresponding to 
25.1 ±1.4 pb -1 , using the same requirements as described above, yielded 20998 W — > ev 
and 1634 Z — > ee candidates; an analysis of 30.7 ± 1.7pb _1 yielded 4516 W — > \iv and 168 
Z — > fi/j, candidates. The spectra are shown in Fig. 2. 
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TABLE II. Analysis results 
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1994-1995 (Preliminary) 
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Nobs 
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Background (%) 


17.3 ± 2.2 


11.0 ±2.4 


17.3 ± 1.1 


10.1 ± 3.7 


Acceptance (%) 


46.1 ± 0.6 


36.3 ±0.4 


22.0 ±0.9 


5.1 ± 0.6 


Efficiency (%) 


66.9 ±4.1 


70.6 ±4.6 


28.6 ± 1.9 


60.9 ± 2.6 


£ (pb" 1 ) 


25.1 ± 1.4 


25.1 ± 1.4 
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30.7 ± 1.7 



The total backgrounds estimated for these event samples are shown in the spectra as 
hashed areas in Fig. 1 and are listed as a percentage of the observed number of events in 
Table I. A major background to the W — > ev sample was from QCD multijet events where 
a jet was misidentified as an electron. It was estimated from data by measuring the Jfif 
distribution of a background-dominated sample, obtained by selecting events containing an 
EM cluster which failed at least one of the electron criteria (isolation, shower shape, and 
track match). We extrapolated this Jfif distribution into the signal region (Jfif > 25 GeV) by 
normalizing the number of events in the background sample to that in the candidate sample 
(without the Jfif requirement imposed) in the region of small Jfif (0 < $t < 15 GeV). The 
backgrounds due to W — > tv — > evvv decay and Z — > e + e~ where one of the electrons 
was lost, were estimated using a Monte Carlo simulation. The multijet background in the 
Z — > ee sample was due to jet-jet or photon-jet events where the jets faked electrons in 
the detector. The amount of this background was estimated by fitting the invariant mass 
spectrum of the Z — > ee events to the sum of the predicted Z boson mass distribution and 
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the experimentally determined multijet background shape. The invariant mass distributions 
lor the jet-jet and photon-jet events were measured separately and then combined to obtain 
the overall multijet background shape. 

The multijet background in the W — > \iv and Z — > fi + \i~ samples was estimated by 
comparing the distribution of energy in the calorimeter between the cones of radii of 0.2 
and 0.6 around the muons with that measured for events containing a non-isolated muon 
and jets. The background in both the W — > \iv and Z — > fi + \i~ samples arising from 
W — ^ tv and Z — > tt decays, as well as the background in the W — > \iv sample arising from 
Z — >■ fi + fi~ where one of the muons was lost were estimated using a Monte Carlo simulation. 
The cosmic ray and random hit backgrounds to the W — > \iv and Z — > fi + fi~ samples were 
estimated from the distributions of muon time of origin relative to the beam crossing. 

Finally, in determining the Z — >■ // cross section, a correction (which is listed as a back- 
ground in Table I) was made for the Drell-Yan process where the lepton pair was produced 
via a virtual photon. This correction was sensitive to the choice of Z mass window. 

The kinematic and geometric acceptances (Tables II) were calculated using a Monte Carlo 
simulation which modeled the detector fiducial volume as well as the measured detector 
resolutions. The calculation used the CTEQ2M [11] pdf and a NLO calculation of p^ and 
Pj by Arnold and Kauffman [5]. The largest contribution to the systematic error in the 
acceptance (Table II) arose from the choice of pdf. We estimated this uncertainty from 
the spread among the values obtained with CTEQ2 [11], MRS [12], and GRV [13]. Other 
errors included were from varying the W mass, the simulation of the p^ and pj distributions, 
radiative corrections, the detector simulation of the Jfif distributions and the detector energy 
scale. The net detection efficiency (Table II) includes both the trigger and offline efficiencies. 
These were estimated from the data using Z — >■ // events since the trigger required only one 
lepton. The electron channel trigger was found to be > 95% efficient; and the muon trigger 
efficiency was 40% (70%) efficient for W(Z) boson events. 

The luminosity was measured by Level trigger scintillator hodoscopes [14] mounted 
at z = ±1.4 m. The north-south coincidence rate was measured and corrected for multiple 
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TABLE III. Cross Section Results for electron (e) and muon (p) channels. When two errors 
are given the first is the statistical error and the second is total systematic error. 

a w ■ B(W ± -» / ± //) (nb) a z ■ B(Z l+l~) (nb) 

1992-1993 

e 2.36 ±0.02 ±0.15 0.218 ± 0.008 ± 0.014 

fi 2.09 ±0.06 ± 0.25 0.178 ± 0.022 ± 0.023 



1994-1995 (Preliminary) 



2.24 ±0.02 ± 0.20 0.226 ± 0.006 ± 0.021 

1.93 ±0.04 ± 0.20 0.159 ±0.014 ±0.022 



Standard Model 2.42±g;^f 0.226^;^ 




PQ 0.20 
D N 0.15 



ev,ee [iv,\i\i ev,ee (av,|j,(j, ev,ee (av,(j,|j, 

FIG. 3. a ■ B for inclusive W and Z boson production. D0 1994-1995 results are preliminary. 
The error bars indicate the combined statistical and systematic errors including the luminosity 
uncertainty. The solid lines are the predicted values calculated using the CTEQ2M pdf and the 
shaded bands indicate the uncertainty in the predictions. 
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interactions. The visible cross section was calculated to be oxo = 46.7 ± 2.5 mb, which 
resulted in a 5.4% relative error on the luminosity determination. This calculation was 
based on an average of the published CDF [15] and E710 [16] measurements of the total, 
elastic, and single diffractive cross sections, with the MBR [17] and Dual Parton Model 
DTUJET-93 [18] Monte Carlo routines used to determine the hodoscope acceptance. 

The resulting cross sections, calculated using Eq. 1, are listed in Table III, where the 
first error given is statistical and the second is the total systematic error, including the 
luminosity uncertainty. These values are compared to the theoretical prediction (taken 
from Ref. [9]) in Fig. 3, together with the CDF results [19]. The total cross sections were 
calculated to be aw = 22.35 nb and oz = 6.708 nb using a numerical calculation program 
from Ref. [1] and using the CTEQ2M pdf [11], M z = 91.19 GeV/c 2 [20], M w = 80.23 ±0.18 
GeV/c 2 [21], and sin 2 9 W = 1 - (M w /M z ) 2 = 0.2259. The branching ratios used are 
B(W -> Iv) = (10.84 ± 0.02)% (calculated following Ref. [22] but with the above M w ), and 
B(Z — > 11) = (3.367 ± 0.006)% [20]. The width of the band in Fig. 3 indicates the error in 
the predicted value, due primarily to the choice of pdf (4.5%) and to the use of a NLO pdf 
with the NLLO calculation (3%) [23]. Figure 4 shows a ■ B for inclusive W and Z boson 
production cross sections from the D0 1992-1993 data, the CDF measurements [19] and the 
measurements at the CERN pp collider [24] as a function of center of mass energy. Good 
agreement between the theoretical prediction and the measurements represents a success of 
perturbative QCD calculations. 

THE TRANSVERSE MOMENTUM DISTRIBUTION 

The large W/Z data samples collected during the 1992-1993 run allow significant im- 
provement in the precision of the da / dpx measurement over previous measurements [25]. 
In this section we describe new, high precision measurements of the px distributions of W 
and Z bosons. The data samples used for the p^ and pj measurements are identical to the 
sample described above for the a ■ B(W — > ev) and a ■ B(Z — > ee) measurements. The p^ 
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FIG. 4. a ■ B for inclusive W and Z boson production as a function of center of mass energy. 
The error bars indicate the combined statistical and systematic errors including the luminosity 
uncertainty. The solid lines indicate the uncertainty of the Standard Model predictions. 
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FIG. 5. Raw distribution (triangles) and multijet background (circles). 

was determined from the hadronic recoil of the W, while the pj was determined from the 
sum of two electron transverse momenta. 

The measurement of the px distribution requires the knowledge of the total amount 
of background, which is listed in Table I, and its shape as a function of the boson pj. 
The shape of the background in the W — > ev sample was obtained by subtracting the p^ 
distribution obtained for a set of very clean electron identification cuts from a px distribution 
of background-dominated sample while accounting for the relative efficiency loss between the 
two cuts. Figure 5 shows the obtained distribution superimposed on the raw p^ distribution. 
The shape of the multijet background in the Z — > ee sample was obtained from data by 
studying the product of the isolation variables of the two electrons as a function of the 
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FIG. 6. Product of isolation variables of the two Z boson electrons versus invariant mass, for 
the standard (top) and loose (bottom) electron identification requirements. 

e + e~ invariant mass, shown in Fig. 6 for the standard and loose electron identification 
requirements. The events from region B (75 < M ee < 105 GeV/c 2 , and Isol • Iso2 > 0.006), 
marked in Fig. 6, were used to parametrize the shape of the multijet background. 

The systematic uncertainties of the measurement arose from uncertainties of: 
i) hadronic energy scale, ii) underlying event contribution, Hi) hadronic resolution, and 
iv) background shape and magnitude. The hadronic energy scale was determined by balanc- 
ing the Z boson px determined from the hadronic recoil and from the transverse momenta 
of the two electrons along the bisector of the angle subtended by the two. The uncertainty 
on the hadronic scale was thus controlled by the number of observed Z candidates and for 
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FIG. 7. Background subtracted pj distribution of data for | rj |< 1.1 (triangles) with smeared 
theoretical prediction [5] (histogram). 

this data sample it produced ~ 20% uncertainty in the measurement. The magnitude of the 
underlying event contribution was also obtained from the Z sample by matching the Z boson 
Pt resolution between data and Monte Carlo and was estimated to be of the order of 10% for 
this measurement. The uncertainty of the background shape and the magnitude was small 
in the low px region but dominant in the high px region. The statistical uncertainty for this 
sample was of the order of 5% per bin (high momentum bins had larger uncertainty, 10 - 
30%). Thus the uncertainty of the pj 7 measurement was dominated by systematic effects, 
most of which were directly controlled by the number of observed Z bosons. 

The statistical uncertainty of the pj measurement was of the order of (10 — 20)% up to 
70 GeV/c. The systematic uncertainties arose from uncertainties of: i) electron energy scale, 
ii) electron energy resolution, and Hi) electron angular (6 and cf>) resolutions. In the present 
analysis, the statistical uncertainty is the dominant uncertainty of the pj measurement. 

Figures 7 and 8 show the comparison between the background subtracted p^ and pj 
distributions, respectively, with theoretical predictions smeared by detector resolutions. The 
data tend to peak at a slightly higher value of px than do the predictions. 
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FIG. 8. Background subtracted pj distribution (solid dots) with smeared theoretical predic- 
tions by [5] (gray histogram) and [6] (dark histogram) superimposed. Top plot: linear scale; Bottom 
plot: logarithmic scale. 
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THE W CHARGE ASYMMETRY 



The W production in pp collisions at yfs = 1.8 TeV is dominantly from a valence- valence 
or valence-sea quark-antiquark interaction. Therefore a VK^VK") is produced primarily by 
the interaction of a u(d) quark from the proton and a d(u) quark from the antiproton. In 
the proton the u valence quark momentum distribution, u(x), is harder than the d valence 
quark distribution, d(x) and, therefore, a VK^VK") is produced with a boost in the proton 
(antiproton) direction. Thus a measurement of the W + and W~ rapidity distributions 
(Yw+-) gives information on parton distribution function (pdf) in the region of low x and 
high g 2 (~ M w ) [26,27]. Because there is a twofold ambiguity in reconstructing Yw in a 
W — > £v decay (due to the fact that the component of neutrino momentum along the beam 
direction is not measured) we measure the Yw distribution indirectly via the charged lepton 
rapidity distribution (Yt), which is a sum of he W rapidity and the lepton rapidity (1^, CM ) in 
the W rest frame: Y#+ = Y\y+ + Y^. M , where Y f GM is determined by the V — A couplings. At 
y/s = 1.8 TeV the asymmetry due to u(x) and d(x) is larger than that from the V — A effect 
and of the opposite sign. The experimentally convenient quantity is the charge asymmetry 
of the lepton pseudorapidity distribution, 

= da(£ + )/dr] - d<r(£-)/dr} 
[1]> ~ d<r(£+)/dr} + da{£-)/df ] ' [ ' 

because it is insensitive to acceptance corrections. Furthermore, because A(—r)) = —A(r)) 
by CP invariance the result can be shown as A(|?y|). A measurement of A(rj) with \rj\ <~ 1.7 
can provide information about the pdfs in the region of x ~ 0.007 — 0.24. 

We present a preliminary result of the W charge asymmetry using W ^ fj,v decays 
observed from the 6.5 pb -1 data of the 1992-1993 run and the first ~ 30 pb -1 data of 
the 1994-1995 run. The data sample was obtained with a single muon trigger: a muon 
with \r)\ < 1.7 and p^ > 15 GeV/c. Additional track quality cuts, identical to those in 
the inclusive W cross section analysis, and pj > 20 GeV/c were imposed offline. For the 
1992-1993 run, 60% of the data were taken with the muon toroid polarity in the forward 
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FIG. 9. D0 Preliminary W decay muon charge asymmetry. The lines correspond to the theo- 
retical predictions using several recent pdfs. 
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direction and the remaining with reversed polarity, while the polarity was flipped every week 
for the 1994—1995 run to minimize possible detector charge asymmetry effects. 

Muon charge misidentification dilutes the charge asymmetry. This systematic effect 
was estimated from the number of same sign pairs in the Z — > fifi sample. The charge 
misidentification probability was 8.8 ± 5.2% for the 1992-1993 run and 2.7 ± 1.5% for the 
1994-1995 run. In addition, if the detector has different acceptance for /i + and it can 
bias the charge asymmetry. Flipping the polarity of the muon toroid compensates this effect. 
The remaining uncompensated luminosity (~ 20% in the 1992-1993 run and ~ 1.5% in the 
1994-1995) was corrected for this bias using the factor derived from the data taken with 
magnet polarities in the forward and reverse directions. 

Figure 9 shows a preliminary W charge asymmetry measurement. Data are compared 
with a leading order calculation with input (y) spectrum obtained from the next-to- 
leading order resummation calculation of Ref. [6]. The data are consistent with the theoret- 
ical predictions with several recent pdfs. 

CONCLUSION 

We have presented the inclusive production cross sections times leptonic branching frac- 
tions in both electron and muon channels. Good agreement between the theoretical calcu- 
lation and our measurements indicates a success of perturbative QCD. We have shown the 
preliminary measurements of transverse momentum distributions of W and Z bosons in the 
electron channel. The large W/Z data samples we obtained improve the precision of the 
da / dpx measurements over previous measurements. A preliminary measurement of the W 
charge asymmetry using W — > \iv decays has also been presented. 
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